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LOCALIZATION OF LAPLACIAN EIGENFUNCTIONS IN 
CIRCULAR, SPHERICAL AND ELLIPTICAL DOMAINS 

B.-T. NGUYEN* AND D. S. GREBENKOV*t t § 

Abstract. We consider Laplacian eigenfunctions in circular, spherical and elliptical domains in 
order to discuss three kinds of high-frequency localization: whispering gallery modes, bouncing ball 
modes, and focusing modes. Although the existence of these modes was known for a class of convex 
domains, the separation of variables for above domains helps to better understand the "mechanism" 
of localization, i.e. how an eigenfunction is getting distributed in a small region of the domain, and 
ff^ decays rapidly outside this region. Using the properties of Bessel and Mathieu functions, we derive 

I the inequalities which imply and clearly illustrate localization. Moreover, we provide an example of 

a non-convex domain (an elliptical annulus) for which the high-frequency localized modes are still 
present. At the same time, we show that there is no localization in most of rectangle-like domains. 
This observation leads us to formulating an open problem of localization in polygonal domains and, 
O more generally, in piecewise smooth convex domains. 

<D 

P^ Key words. Laplacian eigenfunctions. Localization, Bessel and Mathieu functions. Diffusion, 

^ -^ Laplace operator 

II 1. Introduction. A hundred years ago, Lord Rayleigh documented an interest- 

4^ ing acoustical phenomenon that occured in the whispering gahery under the dome of 

^ Saint Paul's Cathedral in London [1] (see also [2l|3]). A whisper of one person prop- 

Ph agated along the curved wall to another person stood near the wall. This acoustical 

i—i effect and many related wave phenomena can be mathematically described by Lapla- 
cian eigenmodes satisfying —Au = A^x in a bounded domain, with an appropriate 

^^ boundary condition: 

^Vl '^ = (Dirichlet), 

1/^ — =0 (iNeumannj, q J^^ 

Cn du 

^ —^hu = (Robin), 
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. . where /i > is a positive constant, and d/dn is the normal derivative directed out- 

^ wards the boundary. It turns out that the eigenmodes that are "responsible" for the 

k^ whispering effect, are mostly distributed near the boundary of the domain and almost 

'^ zero inside. Keller and Rubinow discussed these so-called whispering gallery modes 

Cd and also bouncing ball modes. The existence of such localized eigenmodes in the limit 

of large eigenvalues was shown for any two-dimensional domain with arbitrary smooth 

convex curve as its boundary (so-called high-frequency or high-energy localization) [4 . 

A further semiclassical approximation of Laplacian eigenfunctions in convex domains 

was developed by Lazutkin [5[[6l|7l|8] (see also [9l HOl HT] |12] ) . Chen and co-workers 
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analyzed Mathieu and modified Mathieu functions and reported another type of local- 
ization named focusing modes [13] . These and other localized eigenmodes have been 
intensively studied for various domains, named quantum billiards [T5l fT6| IT7( ITSt [T9] . 
It is also worth mentioning that low- frequency localization of Laplacian eigenfunctions 
in simple and irregular domains has attracted a considerable attention during the last 
two decades [20i i21|i22i i2^ ^ i5j . 

The aim of this paper consists in revisiting and illustrating the aforementioned 
three types of high-frequency localization. For this purpose, we consider circular, 
spherical and elliptical domains for which the separation of variables reduces the anal- 
ysis to the behavior of special functions. Using the properties of Bessel and Mathieu 
functions, we derive the inequalities that clearly show the existence of infinitely many 
localized eigenmodes in circular, spherical and elliptical domains with Dirichlet, Neu- 
mann or Robin boundary condition. More precisely, we call an eigenfunction u of the 
Laplace operator in a bounded domain 1^ C M'^ Lp- localized {p > 1) if it is essentially 
supported by a small subdomain Q^ C 1^, i.e. 

\\u\\Lpin\n^) ^dj^a) ... .-. ^. 

where ||.||lp is the L^-norm, and /i^ is the Lebesgue measure. We stress that this 
"definition" is qualitative as there is no objective criterion for deciding how small these 
ratios have to be. This is the major problem in defining the notion of localization. 
For circular, spherical and elliptical domains, we will show in Sect. [2] and [3] that both 
ratios can be made arbitrarily small. In other words, for any prescribed threshold e, 
there exist a subdomain Via, and infinitely many eigenfunctions for which both ratios 
are smaller than e. Most importantly, we will provide a simple example of a non- 
convex domain for which the high-frequency localization is still present. At the same 
time, we will show in Sect. |4.1| the absence of localization in most of rectangle-like 
domains. This observation will lead us to formulating an open problem of localization 
in polygonal domains and, more generally, in piecewise smooth convex domains. 

2. Localization in circular and spherical domains. 

2.1. Eigenfunctions for circular domains. The rotation symmetry of a disk 
1^ = {x G M^ : \x\ < K\ oi radius R leads to an explicit representation of the 
eigenfunctions in polar coordinates: 

/x^N fcos(n(z?), ^ = 1, ,^,. 

w„fei(r,^) = J„(a„fer/i?)xr ; ^''' . ' , , ^, (2.1) 

I sm(n(/9j, z = 2 (n 7^ 0), 

where Jn[,z) are the Bessel functions of the first kind [26l [27l [28] and ocnk are the 
positive zeros of Jn(z) (Dirichlet), J'^{z) (Neumann) and J!^{z)-\-hJn{z) (Robin). The 
eigenfunctions are enumerated by the triple index nki^ with n = 0, 1, 2, ... counting the 
order of Bessel functions, k = 1, 2, 3, ... counting the positive zeros, and i = 1, 2. Since 
uok2{r,ip) are trivially zero, they are not counted as eigenfunctions. The eigenvalues 
Xnk = <^nfc/^^' which are independent of the last index i, are simple for n = and 
twice degenerate for n > 0. In the latter case, the eigenfunction is any nontrivial 
linear combination of Unki and Unk2- As we will derive the estimates that will be 
independent of the angular coordinate (/?, the last index i will be omitted. 



2.2. Whispering gallery modes. The disk is the simplest shape for ihustrating 



the whispering gahery and focusing modes. The exphcit form (2.1) of eigenfunctions 
ahows one to derive accurate bounds, as shown below. When the index k is fixed, while 
n increases, the Bessel functions Jni^nki^/R) become strongly attenuated near the 
origin (as Jn{z) ~ {z/2)^/n\ at small z) and essentially localized near the boundary, 
yielding whispering gallery modes. In turn, when n is fixed while k increases, the 
Bessel functions rapidly oscillate, the amplitude of oscillations decreasing towards to 
the boundary. In that case, the eigenfunctions are mainly localized at the origin, 
yielding focusing modes. These qualitative arguments are rigorously formulated in 
the following 

Theorem 2.1. Let D = {x e R'^ : |x| < R} be a disk of radius R > 0, 
and Dnk = {x ^ R? : |^| < Rdn/oink}, where d^ = n — v?!'^ , and ank ctre the 
positive zeros of Jn{z) (Dirichlet) , J!^{z) (Neumann) or J!^{z) + hJn{z) for some 
h > (Robin), with n = 0,1,2,... denoting the order of Bessel function Jn{z) and 
k = 1,2,3, ... counting zeros. Then for any p > I (including p = oo^^ there exists a 
universal constant Cp > such that for any k = 1,2,3,... and any large enough n, 
the Laplacian eigenfunction Unk for Dirichlet, Neumann or Robin boundary condition 
satisfies 

'^^"^-^^-^ <C,nH^2-^^V3. (2.2) 

\\Unk\\Lp{D) 

This estimate implies that 

j.^ WUnkUiP..) ^ Q^ ^j^.j^ j.^ ^2{Dnk) ^ ^ (2.3) 

The theorem shows the existence of infinitely many Laplacian eigenmodes which are 
I/p-localized near the boundary dD (see Appendix [A| for a proof). In fact, for any 



prescribed thresholds for both ratios in (1.2), there exists no such that for all n > no 



the eigenfunctions Unk are L^-localized. These eigenfunctions are called "whispering 



gallery eigenmodes" and illustrated on Fig. 2.1 

A simple consequence of the above theorem is 

Corollary 2.2. For any p > 1 and any open subset V compactly included in D 
(i.e., VndD = $), one has 



hm f!^^^ = 0. (2.4) 



As a consequence. 



C,(F).infm^^ (2.5) 

In fact, for any open subset V compactly included in D, there exists no such that for 



ah n > no, F C D^k so that Wunkh^iv) < Wunkh^iDr^k) yielding Eq. ( |2.4 ). 

In the same way, the localization also happens for any circular sector. 

In three dimensions, the existence of whispering gallery modes in a ball follows 
from the following 

Theorem 2.3. Let B = {x eR^ : \x\ < R} be a ball of radius R, and Bnk = 
{x ^ B : < |x| < Rsn/oLnk}, where s^ = (n + 1/2) — (n + 1/2)^/^ and ank ^^^ 
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Fig. 2.1. Formation of whispering gallery modes u^k in the unit disk with Dirichlet boundary 
condition: for a fixed k (k = for top figures and k = 1 for bottom figures), an increase of the index 
n leads to stronger localization of the eigenfunction near the boundary. 



the positive zeros of jn{z) (Dirichlet), j^(^) (Neumann) or j'^{z) + hjn{z) for some 
h > (Robin), with n = 0, 1, 2, ... denoting the order of the spherical Bessel function 
jn{z) o^nd k = 1,2,3,... counting zeros. Then, for any p > 1 (including p = oo^^ 
there exists a universal constant Cp > such that for any k = 1,2,3,... and any 
large enough n, the Laplacian eigenfunction Unk 'with Dirichler, Neumann or Robin 
boundary condition satisfies 



J^nk\\Lp{Br,k) 



\\Unk\\ 



LpiB) 



<Cp(n + l/2) 



i + _ 



exp 



1/3N 



(n>l). (2.6) 



^45 a consequence. 



lim 



\Unk\\Lp{Br,k) 
WUnkhpiB) 



0, 



while lim 



fJ^si^nk) 
M3(B) 



1. 



(2.7) 



As for the disk, the above results show that infinitely many high-frequency eigen- 
functions are Lp-localized near the boundary of the ball (see Appendix [b| for a proof). 

2.3. Focusing modes. The localization of focusing modes at the origin is de- 
scribed by 

Theorem 2.4. For each R e (0, 1), let D{R) = {x eR'^ : R < \x\ < 1}, and D 
be the unit disk. Then, for any n = 0,1,2,...^ the Laplacian eigenfunction Unk w^th 
Dirichlet, Neumann or Robin boundary condition satisfies 



lim 



\\Unk\\ 



Lp(DiR)) 



(1 



\\Unk\\ 



Lp(D) 



j^2-p/2y/p 





(l<p<4), 
(p>4). 



(2.^ 



The theorem states that for each non-negative integer n, when the index k increases, 
the eigenfunctions Unk become more and more L^-localized near the origin when p > 4 
(see Appendix [A| for a proof). These eigenfunctions are called "focusing eigenmodes" 
and illustrated on Fig. |2.2[ The theorem shows that the definition of localization 
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Fig. 2.2. Formation of focusing modes u^k ^^ the unit disk with Dirichlet boundary condition: 
for a fixed n (n = for top figures and n = 1 for bottom figures)^ an increase of the index k leads 
to stronger localization of the eigenfunction at the origin. 



is sensitive to the norm: the focusing modes are L^-locahzed for p > 4 (including 
p = oo), but they are not L^-locahzed for p < 4. In fact, as the amphtude of 
oscihations of the focusing modes exhibits a power law decay from the origin towards 
the boundary (see Fig. 2.2), the p values control the behavior of the L^-norm and 
determine whether the ratio of these norms vanishes or not in the limit k ^ oo. 

A similar theorem can be reformulated for a ball in three dimensions. 

Theorem 2.5. For each R e (0, 1); let B{R) = {x eR^ : R < \x\ < 1}, and 
B be the unit ball. Then, for any n = 0,l,2,...^ the Laplacian eigenfunction Unk with 
Dirichlet, Neumann or Robin boundary condition satisfies 



k^oo \\Unk\\Lp{B) 



{1-R^-py/p (1 <p<3), 

(p>3). 



(2.9) 



(see Appendix [B] for a proof). 

3. Localization in elliptical domains. 

3.1. Eigenfunctions for elliptical domains. It is convenient to introduce the 
elliptic coordinates as 



xi = acoshrcos6>, 
X2 = asinhrsin^. 



(3.1) 



where a > is the prescribed distance between the origin and the foci, r > and 
< ^ < 27r are the radial and angular coordinates. A filled ellipse is a domain 
with r < R so that its points (xi,X2) satisfy xf/A^ + x^/B'^ < 1, where R is the 
elliptic radius and A = a cosh i? and B = asinhi? are the major and minor semi- 
axes. In the elliptic coordinates, the separation of the angular and radial variables 
leads to Mathieu and modified Mathieu equations, respectively [13l |14- Periodic 
solutions of the Mathieu equation are possible for specific characteristic values c. 
They are denoted as cen{0,q) and sen-\-i{0^q) (with n = 0,1,2,...) and called the 
angular Mathieu functions of the first and second kind. Each function ce^(^,g) and 
sen+i(^,^) corresponds to its own characteristic value c (the relation being implicit, 
see [in]). 
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For the radial part, there are two Unearly independent solutions for each char- 
acteristic value c: two modified Mathieu functions Mc^^-' (r^q) and Mc}^ ^ (r, g) corre- 
spond to the same c as cen(6^, ^), and two modified Mathieu functions MsJ^^;^(r, q) and 

MsJ^^]^(r, q) correspond to the same c as sen+i(^, q). As a consequence, there are four 
families of eigenfunctions (distinguished by the index i = 1, 2, 3, 4) in a filled ellipse: 

Unki = ceniO,qnki)Mcll\r,qnki), 

Unk2 = Cen{0,qnk2)Mcl^\r,qnk2), 
UnkS = Sen+l(6>,gn/c3)Ms^^|i(r,g'nfc3), 

Unki = se„+i(6',g„fe4)MsJj|i(r,(2„fc4), 



(3.2) 



where the parameters Qnki are determined by the boundary condition. For instance, 
for a filled ellipse of radius R with Dirichlet boundary condition, there are four indi- 
vidual equations for the parameter q„fei, for each n = 0,l,2,...: 



McW(i?,g„fei) = 0, Mc(,2)(i?,5„fe2) = 0, 



.(1) 



(2) 



Ms^Yi (^. 9"fe3) = 0, Ms':i^ (^. Inki) = 0, 



(3.3) 



each of them having infinitely many positive solutions qnki enumerated by /c = 1,2,... 
]. The associated eigenvalues Xnki are determined as 



'^nki 



^Qriki 



(3.4) 



The above analysis can be applied almost directly to an elliptical annulus 1], 
i.e. a domain between an inner ellipse Fi and an outer ellipse F2, with the same 
foci. In elliptic coordinates, ft can be defined by two inequalities: Ri < r < R2 and 
< ^ < 27r, where the prescribed radii Ri and R2 determine Fi and F2, respectively. 

We consider two families of eigenfunctions in Q: 



^n/ci — ce^ 



n{0,qnki) ankiMcl^\r,qnki) ^bnkiMc^^\r,qnki) 



Unk2 = Sen+l(6>,gnfe2) «n/c2Ms[^^ji (r, ^^^2) + Kk2Ms^nlli^^ Qnk2) 



(3.5) 



The parameters a^/ci, ^n/ci and qnki (^ = 1,2) are set by boundary conditions and 
the normalization of eigenfunctions. For Dirichlet boudary condition, one solves the 
following equations: 



(3.6) 



Mc^^\Rl,qnkl)Mc'^\R2,qnkl)-Mc'^\R2,q„kl)Mc'^\Rl,q„kl)=0, 
Msi'jl(i?l,g„fe2)Msi'ji(i?2,9nfe2) - Msi'ji(i?2,9„fc2)Msi'ji(i?i,9„fc2) = 0. 

For n = 0, 1, 2, ..., each of these equations has infinitely many solutions qnki enumer 
ated by fc = 1, 2, 3, ... [29. The eigenvalues are determined by Eq. (3.4). 



3.2. Bouncing ball modes. For each a G (O, f), we consider the ehiptical 
sector Vtoc inside an ehiptical domain Vt: 

^cx = {Ri < r < i?2, 6> G (a, TT - a) U (tt + a, 27r - a)} . 



Theorem 3.1. Let ft be a filled ellipse or an elliptical annulus (with a focal 
distance a > 0). For any a e (O,^), p > 1 and i = 1,2,3,4 (for filled ellipse) or 
i = 1, 2 (for elliptical annulus), there exists A^^ > such that for any XnUi > ^a 



\\^nki\\Lp{n\n^) ( 16a \^'^ 



\\^nki\\Lp{n) 

where 



< Dn ( —^ — 7^ j exp i-a\/\nki sin (^- + - j - sina j , (3.7) 




(see Appendix |d] for a proof; a similar exponential bound for the I/2-norm was 
recently derived in [3O4). Given that X^ki ^ 00 as /c increases (for any fixed n and z), 
while the area of flex can be made arbitrarily small by sending a -^ 7r/2, the theorem 
implies that there are infinitely many eigenfunctions Unki which are Lp-localized in 
the elliptical sector Q^- 

lim """^^"^r^"\"°^ = 0. (3.9) 

k^oo WUnkillLpin) 



These eigenfunctions are called "bouncing ball modes" and illustrated on Fig. 3.1 We 
note that similar results were already known for a filled ellipse and, more generally, 
for convex planar domains with smooth boundary [^ [13] . Although our estimates are 
specific to elliptical shapes, they are explicit, simpler and also applicable to Lp norms 
and to elliptical annuli, i.e. non-convex domains. 



The quality of the above estimates was checked numerically. Figure 3.2 shows the 



ratio 1^^^^ L2(n\ncx) ^^^ -^g upper bound for two families of eigenfunctions in a filled 
ellipse and an elliptical annulus. One can clearly see the rapid exponential decay of 
this ratio when k increases that implies the localization in a thin sector around the 
vertical (minor) axis. Note that the upper bound is not sharp and can be further 
improved. 

4. Discussion. The explicit estimates from previous sections provide us with 
simple examples of domains for which there are infinitely many Lp-localized eigen- 



functions, according to the definition (1.2). Most importantly, the high-frequency 



localization may occur in both convex and non-convex domains. This observation 
relaxes, at least for elliptical domains, the condition of convexity that was significant 
for the construction of whispering gallery and bouncing ball modes by Keller and 
Rubinow [4 and for semiclassical approximations by Lazutkin [II |8]. At the same 
time, these approximations suggest the existence of L^-localized eigenfunctions for a 
large class of domains. How large is this class? What are the relevant conditions on 
the domain? To our knowledge, these questions are open. In order to highlight the 
relevance of these questions, it is instructive to give an example of domains for which 
there is no localization. 
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Fig. 3.1. Formation of bouncing ball modes Unki in a filled ellipse of radius R = 1 (top) and 
an elliptical annulus of radii 0.5 and 1 (bottom), with the focal distance a = 1 and Dirichlet boundary 
condition. For fixed n = 1 and i = 1, an increase of the index k leads to stronger localization of the 
eigenfunction near the vertical semi-axis (Kmax = 200, see Appendix\Cy. 



numerical 

bound 




-numerical 
-bound 



numerical 

bound 




-numerical 
-bound 



a^^ 



Fig. 3.2. The ratio 



nkl\\L2(n\na) 



aV 



(solid blue line) and its upper bound 



, , ^^_. ,3.7) (dashed red 

\\'^nkl\\L2(^a) / r- r- ^ k 

line) in a filled ellipse of radius R = 1 and focal distance a = 1 (top) and in an elliptical annulus 
of radii R = 0.5 and R = 1 and focal distance a = 1 (bottom), with n = and a = 7r/4 (left) and 
n = 1 and a = 7r/3 (right). For numerical computation of Mathieu functions, we used Kmax = 200 
(see AppendixlCy. 



4.1. Rectangle-like domains. Rectangle-like domains, 1] = (0,-^i)x...x(0,-^(i) C 
R^ (with the sizes £i > 0), may seem to present the simplest shape for studying the 
Laplacian eigenfunctions as they are factored and expressed through sines (Dirichlet), 
cosines (Neumann) or their combination (Robin): 



(Xi, ...,Xd) = U^Ji^l) ' • • ^nji^d), A, 



id)( 
rid 



XW + ... + X(f), (4.1) 
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with the multiple index ni-.-n^^, and Un-{xi) and An, {i = l,...,(i) corresponding to 
the one-dimensional problem on the interval (0,^i): 

uli\x) = sm{7rnx/£i), X^^ = Tr^nV^- , n = 1,2,3, ... (Dirichlet), 

u^\x) = cos{7rnx / ii) ^ X^-' = Tr^n^/-^^, n = 0, 1,2, ... (Neumann) 

(Robin boundary condition will not be considered here). The situation is indeed 
elementary for rectangle-like domains for which all eigenvalues are simple. 

Theorem 4.1. Let 1] = (0,fi) x ... x (0,-^^^) cR^ be a rectangle-like domain with 
sizes ii > 0, ..., id > such that 

fje'jiQ ^i + j. (4.2) 

f(Q denoting the set of rational numbers). Then for any p > I and any open subset 

vcn, 

C,{V)= inf /|l!^^^^k(n\>o. (4.3) 

ni,...,nd [ \\Um,...,nd\\LpiQ) J 

The proof is elementary (see Appendix IeI) and relies the fact that all the eigenvalues 



are simple due to the condition (4.2). The fact that Cp{V) > for any open subset V 
means that there is no eigenfunction that could fully "avoid" any location inside the 
domain, i.e., there is no L^-localized eigenfunction. Since the set of rational numbers 



has zero Lebesgue measure, the condition (4.2) is fulfilled almost surely, if one would 



choose a rectangle-like domain randomly. In other words, for most rectangle-like 
domains, there is no L^-localized eigenfunction. 

When at least one ratio i'f/i'j is rational, certain eigenvalues are degenerate, and 
the associate eigenfunctions become linear combinations of products of sines or cosines. 
For instance, for the square with ii = £2 = tt and Dirichlet boundary condition, the 
eigenvalue Ai^2 = 1^ + 2^ is twice degenerate, and 1^1,2(^1,^2) = ci sin(a::i) sin(2x2) + 
C2 sin(2xi) sin(x2), with arbitrary constants Ci and C2 {cf -\- c^ ^ 0). Although the 
computation is still elementary for each eigenfunction, it is unknown whether the 



infimum Cp{V) from Eq. (4.3) is strictly positive or not, for arbitrary rectangle-like 
domain Q and any open subset V. The most general known result for a rectangle 
n = {0,£i)x (0, £2) states that C2{V) > for any V C 17 of the form V = {0,£i)x uj, 
where u is any open subset of (0,^2) 131 . Even for the unit square, the statement 
Cp{V) > for any open subset V seems to be an open problem. More generally, 
one may wonder whether Cp{V) is strictly positive or not for any open subset V in 
polygonal convex domains or in piecewise smooth convex domains. To our knowledge, 
these questions are open. 

5. Conclusion. We revived the classical problem of high-frequency localization 
of Laplacian eigenfunctions. For circular, spherical and elliptical domains, we derived 
the inequalities for L^-norms of the Laplacian eigenfunctions that clearly illustrate 
the emergence of whispering gallery, bouncing ball and focusing eigenmodes. We gave 
an alternative proof for the existence of bouncing ball modes in elliptical domains. 
This proof relies on the properties of Mathieu functions and is as well applicable to 
elliptical annuli. As a consequence, bouncing ball modes also exist in non-convex do- 
mains. At the same time, we showed that there is no localization in most rectangle-like 
domains that led us to formulating the problem of how to characterize the class of do- 
mains admitting high-frequency localization. In particular, the roles of convexity and 
smoothness have to be further investigated. The problem of localization in polygonal 
convex domains or, more generally, in piecewise smooth convex domains are open. 
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Appendix A. Proofs for a disk. 



The proof of Theorem 2.1 is based on several estimates for Bessel functions and 
their roots that we recall in the following lemmas. In this Appendix, jj^^k and jl ^ 
denote all positive zeros (enumerated by k = 1,2,3, ... in an increasing order) of the 
Bessel function Ji^{x) and its derivative J^(x), respectively. 

Lemma A.l. For any n = 1,2,3,... and any e G (0,2/3)^ the Bessel function 
Jn{x) satisfies 



< Jn{nz) < 2-^'/^ V z G (0, 1 - n'-i). (A.l) 

Lemma A. 2. The first zeros jn,i and j'^ ^ with n = 1, 2, ... satisfy 



^ < fn,i < Jn,i < Vn + 1 (Vn + 2 + 1) . (A.2) 

Lemma A. 3. For large enough n, the asymptotic relations hold JWj : 



Un) = CW'" + 0(n-^l^) [C', = ^1^ ^ 0.4473) (A.3) 

j;^i = n + n^/^ C^ + 0{n-^^^) (C^ = 0.808618...). (A.4) 

^4^ a consequence, taking smaller constants (e.g., Ci = 0.447 and C2 = 0.8086^)^ one 
gets lower bounds for large enough n: 

Jn(n)>Cin-^/^ (n>l) (A.5) 

j;,i >n + C2n^/^ (n>l), (A.6) 

Lemma A. 4. For fixed k and large v, the Olver^s expansion holds J3^\3^\3^ 

r 1/3 3,2 -1/3 5-(5| _i 479(^^ + 204 -5/3 

•^"^'^ ^ 10 ^ 350 63000 /^ j^ 

20231SI - 27550SI _^i^ ^, _^^ 

8085000 ^ ^' 

where 6k = —ak2~^^^ > and a^ are the negative zeros of the Airy function (e.g., 
61 = 1. 855757. ..y). Taking Ck = Sk -\- e (e.g., e = 1), one gets the upper bounds for j^^k 
for V large enough 

3v,k<y^CkV^''^ (i^>l). (A.8) 



Lemma A. 5. For fixed v and large k, the McMahon^s expansion holds f2Si/ (p. 
506) 

^•^'^ = ^^^>- V2) - 8(,./:(l"_\/2)/2) + ^(Vfc^). (A.9) 
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Lemma A. 6. The absolute extrema of any Bessel function Jy{z) progressively 
decrease 126] (p. 4^8), i.e. 



\U3'u,i)\ > k.(i:,2)i > \mjI,3)\ > - 



(A.IO) 



Lemma A. 7. The k-th positive zero ank of the function J!^{z) + hJn{z) for any 
h > lies between the k-th positive zeros jn,k ci'^d j'^ ^ of the Bessel function Jn{z) 
and its derivative J'^^z): 



Jn,k < ^nk < jn,k' 



(A.ll) 



Proof This is a direct consequence of the minimax principle that ensures the 
monotonous increase of eigenvalues with the parameter h [37 . D 



Using these lemmas, we prove Theorem 2.1 



Proof The proof formalizes the idea that the eigenfunction Unk is small in the 
large subdomain Dnk = {x ^ D : \x\ < Rdn/oink} (with dn = n — n^/^) and large in 
the small subdomain Ank = {x G D : Rn/ank < |^| < Rj'n i/oLnk}- Since A^k C D^ 
we have for 1 < p < oo 

Rdn/ank dn 

lu, jp 117/ Jl^ / dr r\Jn{rank/R)\^ J dzz \Jn{z)\P 

II wn ^11 iiT) 



p ^ J dr r \Jn{rank/R)\P J dzz \Jn{z) 

Rn/ank ^ 



The numerator can be bounded by the inequality (A.l) with e = 1/3: 

J dzz \MzW < (2-^^^/^)' f < 2--^^^/3 !^ (n = 1,2,3,...). 



In order to bound the denominator, we use the inequalities (A.5, A. 6) and the fact 



that Jn{z) increases on the interval [n, j^ i] (up to the first maximum at j^ ;l)* 
J dzz \Uz)\P > |Jn(n)r^:^^^ > [C,n-V3]. {n ^ Cn^p' - 



> Cf C2n(^-^)/^ 



for n large enough, from which 



hnkh^iDnk) ^ n'^^^^P 



+ ^2-^''''/3 



IKkhAAnk) Ci(2C2)Vp 



(n>l) 



that implies Eq. (2.2). 



For p = 00, one has 



II 11 II II max \Jn(z)\ 

WUnkUo^jDnk) W^nkU^jDnk) _ 0<z<dn 



UnkllL^iD) W'^nkh^iAnk) ^^^ \Jn{z)\' 
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Using the same bounds as above, one gets 



max \Jn{z)\ < 2 "^^ ^^^, max \Jn{z)\ > Jn{n) > Cin ^/^ 

0<Z<dn '^<^<fn,l 



that imphes Eq. (2.2). 



Finahy, from Lemmas |A.4| and |A.7[ we have 

2 



1 > 



dl (n - n2/3)2 






(n>l) 



so that the ratio of the areas tends to 1 as n goes to infinity. D 



Let us prove Theorem 2.4 



Proof. For p = oo, the exphcit representation ( |2.1[ ) of eigenfunctions leads to 

(A.12) 



II II max \Jn{(^nkr)\ max \Jn{(^nkr) 



hnkh^iD) max\Jn{ankr)\ \Jn{fn,l)\ 

rG[0,lJ ' 



where we used the fact that the first maximum (at j^ i) is the largest (Lemma A.6). 
Since lim ank = co, the Bessel function Jn{o^nk^) with k ^ 1 can be approximated 

in the interval [R^ 1] as [28J 

/ nir 7r\ 



Jn{(^nkr) 



TTOinkr 



(A.13) 



It means that there exists a positive integer Kq and a constant Aq > (e.g., Aq = S/tt) 
such that 



|J„Kfer)| <^j^^< J^^, Vr e [R,l], k > Ko. 



(A.14) 



Given that the denominator in Eq. (A.12) is fixed, while the numerator decays as 

— 1/2 I I 

a^^ , one gets Eq. (2.8) for p = oo. 
For p > 4, the ratio of Lp norms is 



\\^rik\\L^(^D{R)) ar^kR 



j dr r \Jn{r 



hnkWl^^D) Ydrr \Jn{r)\P 



(A.15) 







The inequality (A.14) allows one to bound the numerator as 

OCnk 



I 



dr r \Jn{r)\P < A 



Otnk 



■P/-2 _ 



^g/2^^2-p/2 _ 1] 

p/2-2 



a. 



2-P/2 

nk ' 



OCnkR OCnkR 

while the denominator can be simply bounded from below by a constant 

Otnk 1 



f dr r \Jn{r)\P > f dr r \J„{r)\P. 
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As a consequence, the ratio of Lp norms in Eq. (2.8) goes to as A: increases. 
For 1 < p < 4, one can write 



hnkfr 



LpiDiR)) 



hnkWl 



1 



LAD) 



fp,n{(^nkR) 
Jp,n\^nk) 



where 



fp,„iz)= drr\Jn{r)\P. 



(A.16) 



As discussed in Remark A.l, the function /p,n(^) behaves asymptoticahy as z 



2-P/2 

for large z^ i.e., there exists < Cp^n < oo such that for any 6: > 0, there exists zq > 
such that for any z > zq [cf. Eq. (A. 17)] 



\^p,r, 



- e)z^-P/^ < U,n{z) < {cp,n + e)z^-P/\ 



from which one immediately deduces 



_^2-p/2 < 



Jp,n \(^nk ) 



< 



:^2-p/2_ 



As a consequence, for any i? < 1, one can always choose s such that the right-hand 
side is strictly smaller than 1 so that the ratio of Lp norms is then strictly posit ive. 
Moreover, the limiting value is 1 — R^~^^'^ that completes the proof of Eq. (2.8) for 
1 <p<4. D 

Remark A.l. For 1 < p < A, the function fp^n{z) defined by Eq. (A.16) 
asymptotically behaves as z^~^l^ for large z, i.e., the limit 



lim 



Jp,n\^) 



(A.17) 



exists, is finite and strictly positive: < Cp^n < oo. Although this result is naturally 
expected from the asymptotic behavior (A. 13) of Bessel functions, its rigorous proof 



is beyond the scope of the paper. An upper bound for the limit (i.e., Cp^^i < co) can be 



easily deduced from the inequality (A. 14). A lower strictly positive bound (i.e., Cp^n > 



0) would require more careful estimations. The most difficult part consists in proving 
the existence of the limit, as the numerical computation of the function /p,n(^)/^^~^^^ 
at large z shows its oscillatory behavior with a slowly decaying amplitude. Since the 
statement of Theorem 2.4 for 1 < p < 4 relies on this conjectural result, Eq. (2.8) 



was also checked numerically and presented on Fig. |A.1[ 

Appendix B. Proofs for a ball. 

In this Appendix, we generalize the previous estimates to a ball B = {x e 
R^ : 1^1 < R} of radius R > 0. We recall that the Laplacian eigenfunctions in 
spherical coordinates are 



Unkl (r, 0,(p) = jn {ankr/R)Pn (cOS 6) 6^^^ , 



(B.l) 



where jn(^) are the spherical Bessel functions of the first kind (not to be confused 
with zeros jn,fe), 



jn{z) 



Jn+l/2{z)^ 



(B.2) 
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Fig. A.l. The ratio —r. ^ as a function ofpforn = 1 and R = 0.8, in two dimensions 

(left) and three dimensions (right). Three color curves correspond to three eigenf unctions with 
k = 100, k = 1000 and k = 10000, whil e the solid black curve shows the theoretical limit as k ^^ oo 
given by Eq. [2^ in 2D and Eq. \2.9]j in 3D. Note a slower convergence (deviations) for p r^ 4 in 
2D andpr^S in 3D. 



Pn{z) are associated Legendre polynomials, and ank are the positive zeros oi jn{z) 
(Dirichlet), Jn{z) (Neumann) or j'j^{z) + hjn{z) for some h > (Robin), which are 
enumerated by k = 1, 2, 3, ... for each n = 0,l,2,.... We will derive the estimates that 
do not depend on the angular coordinates 6 and (/?, so that the last index / will be 
omitted. 

We start by recalling and extending several classical estimates. 

Lemma B.l. For any v G IR+ and any x G (0, 1)^ the Kapteyn^s inequality holds 

IMI 



< Jjy{iyx) < 






(B.3) 



Now we can prove the following 

Lemma B.2. For any v > 1 and < e < 2/3^ one has 



/TT /^2 1 ^ 

< 7,. ^(x) < \ —- exp [ - — -y 



VxG (0,z/-z/'+i/^). (B.4) 



Proof. Using the Kapteyn's inequality (B.3) and taking x = vz with 2: G (0, 1 
^ye-2/3^^ one has 



3y-l/2{^z) 



Vj,{vz) A^„_ .,, „_ z--V2eKi-')''' 



^ <\—f(z), with f{z)- ^ ^ ^ 



Substituting u = Vl — ^^ ^ (0? 1)? o^^ g^^^ 

u/2 



f{z) = 



(1 



,2\^-277 ^211 



{i + uY 



62 



(i^uY 



l-U 



'2u 



-1/4 



Using the inequality 



1 + ii 3 
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one gets 



/w< 



1- -w" 



1-1/4 



< 62 



1 U^ 

3 



(1 + 0)^ 
Since z < 1 — z/^~^/^, one has 

u = \/l - z'^ > \J\- z > iy^~i > ui~i , 
from which 



1-1/4 



/W < e^ 



l--i/^ 



-1/4 



< 62 



1- -U 



e-1 



Since 



(l-x)i<6"\ VxG(0,1), and < -z/'"^ < - < 1, 



one finally gets 



/(z)<exp(i ^ ^/^^ --^-^^ ' ^--1 



1 1 



2 tz/i-^ J ^V26 3/ ^V33 



that completes the proof. D 

As a consequence, taking v = n -\- 1/2 and e = 1/3, one has 

Lemma B.3. For n = 1,2,... and any z e (0,n + 1/2 - (n + 1/2)^/^), 



jn{z) < 



7T / 2 1 

exp - — - \ n -\- - 

2n + 1 ^ I 3 3 V 2 



l/3^ 



(B.5) 



The lemmas for Bessel functions and their zeros from Appendix [A] allow one to 
get similar estimates for spherical Bessel functions jn(^), their positive zeros jri,k and 
the positive zeros 7^ ^ of j^(z). They are summarized in the following 

Lemma B.4. For n large enough, 

(B.6) 
(B.7) 
(B.8) 
(B.9) 



jn{n + 1/2) >(7i(n + 1/2) "^/e (Cj = v/;72 Ci), 
7„,fc < (n + 1/2) + 4(n + 1/2)1/3^ 
7;_i > n + 1/2 + C2(n + 1/2)^/3 (Ca = 0.80), 

7n,fc < Oink < ln,k- 



Proof. Prom Lemma [A. 3[ we have 

jv-l/2{v) -- 



2t^'^''^''^^yi^"'~'^' = ^"'~'^'' 



from which (B.6) follows by taking u = n -\- 1/2. 



The zeros 7^,^ of the spherical Bessel function jn{z) are also the zeros of the Bessel 
function Jn+i/2(^) so that (B.7) follows directly from Eq. (A. 8) for z^ = n + 1/2 large 
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enough. 

The inequahty (B.8) fohows from the asymptotic expansion of 7^ ^ for large n [28] (p. 

441) 

7;^! = n + 1/2 + 0.8086165(n + 1/2)^/^ - 0.236680(n + 1/2)-^/^ 
- 0.20736(n + 1/2)-^ + 0.0233(n + 1/2)-^/^ + ... (n > 1). 



Taking C2 = 0.80, one gets the inequahty (B.8). 

Finahy, the inequahties (B.9) fohow from the general minimax principle as for the 

disk. D 

We also prove that the first maximum of the spherical Bessel function at 7^ i is 
the largest (although this is a classical fact, we did not find an explicit reference). 



Lemma B.5. For an integer n > 0, one has 

max 7n(^) = jn(7n,l)- 
a?G(0,oo) 



(B.IO) 



Proof. The spherical Bessel function jn{x) satifies 

x^f; + 2xj; + [x^ - (n + l)n]jn = 0. 
Denoting k. = n{n + 1), one can rewrite this equation as 

Jr^yX) — 2 •) 



from which 



dx 



^ifni^)f 



d_ 
dx 



ifni^)) 



A- {fnix)f 



2{x^-K)jUx)j':{x)-2x{j'^{x)f 



'^fnix) 



{x^ — k) 
[2xj'„{x)+ (a;^ - k) jn{x)] 



2x2 



-X{x)f:ix) 



2i;(a;)j;'(x) 

2xhi 



2x 



{jUx)f[2{x'-K) + K] 



2XK: 

2x 



d 



Un{^)f 



{jUx)r [2x' - k] 



dx ^^"(")] • 



Now, if we put 



then 



K{x) = Jn(^) + 



dx 



An{x) = 



2x 



(x^ — K,y 



^Afni^)f 



{j'Jx)f[2x'-n\<0 



for all X > xi ^ ^ i-^- ^n(^) monotonously decreases. Given that An{l'n^k) 



jUln,k) and 



n{n + 1) 



< 7n,l < 7n,2 < 
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we get the conclusion. D 

Now, we can prove Theorem |2.3[ 

Proof. As earher, the proof formahzes the idea that the eigenfunction Unk is smah 
in the large subdomain Bnk = {x G B : \x\ < Rsn/o^nk} (with Sn = {n-^ 1/2) — (n + 
1/2)^/^) and large in the small subdomain Ank = {x ^ B : R{n -\- l/2)/ank < \x\ < 
^7n i/(^nk}' Since Ank C 5, we have for 1 < p < oo 

RSn/ank Sn 

lu, jip IL, jip / rfr r2 |j„(ra„fe/i?)|P J dzz"^ \jn{z)\P 



l|Wnfc||>, fR) \\u„k\\r (A tl «7;,i/""fc T«,i 

i?(n+l/2)/a^fc n+1/2 



The numerator can be bounded by the inequality (B.5): 



^ ' ' (.1 = 1,2,3,...). 



In order to bound the denominator, we use the inequalities (B.6, B.8) and the fact 
that jn{z) increases on the interval [n + 1/2,7^ ;l] (^P ^o ^^^ ^^^t maximum at 7^^i): 

/ (i2;2; |jn(^)r > bn(n + 1/2)]^ '- ^ 

n+1/2 

> [(7,(n + 1/2)-/^]^ (n + 1/2 + (7.(n+l/2)V3)3_(n+ 1/2)3 ^ ^^^^^^ ^ ^/2)V3-5./e 
for n large enough, from which 



< ~ ~ ^exp 



{l-l{^+\) '^') (n+l/2)i/3+^/(3^) (n » 1) 



lknfe||L,(A„,) Ci{3C2y/P 

that implies Eq. (|2.6|). The case p = oo is treated similarly. Finally, from Lemma 



B.4 we have for n large enough 



M3(gnfc) ^ / M . A > (n + 1/2 -(n+ 1/2)^/3)3 
^is{B) \a^J ^^3^^ ^ (n+1/2 + 4(n + 1/2)1/3)3 

SO that the ratio of volumes tends to 1 as n goes to infinity. D 



The proof of Theorem |2.5| for a ball is similar to that of Theorem 2.4 



Proof. For p = 00, the explicit representation (B.l) of eigenfunctions leads to 

11 ,, max \jn{(^nkr)\ max \jn{(^nkr)\ 



hnkh^iB) max \jn{ankr)\ Jn(7n,l) 
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where we used the fact that the first maximum (at 7^ i) is the largest (Lemma B.5). 
Since lim a^k = 00, the spherical Bessel function jn{<^nk^) with k ^ 1 can be 

approximated in the interval [R, 1] as [28] 



jn{0Lnkr) = 



'^(^nkr 



Jn+l/2{(^nkr) 



1 / n + lTT 

COS ankT 



(B.12) 



It means that there exists a positive integer Kq and a constant Aq > (e.g., Aq = 2) 
such that 



\jn{ankr)\ < ^^ < ^^^ yre[RA], k>Ko. 

C^nkr OLnkR 



(B.13) 



Given that the denominator in Eq. (B.ll) is fixed, while the numerator decays as 
a~l, one gets Eq. (2.9). 

For p > 3, the ratio of Lp norms is 



ll^^fe|lLp(D(i?)) _ ar^kR 

11?/ / F ~ "-'^ 





The inequality (B.13) allows one to bound the numerator as 

ank 



ank 



dr r^ \jn{r)\P < Al I dz z^-p = ^^^^—^ a^'^ 



p-3 



nk ' 



OinkR OinkR 

while the denominator can be simply bounded from below by a constant 

<^nk 1 

drr^ \jn{rW>Jdrr^ \Ur)\P. 



As a consequence, the ratio of Lp norms goes to as /c increases. 
For 1 < p < 3, one can write 



ll^^fe|lLp(B(i^)) _ _ fp^njc^nkR) 

II II r) -L ~ , . • 



11'^^^ 1 1 Lp(S) 



Jp,n\^nk) 



where 



z 

U,^{z)^ Jdrr^ |j„(r)|^ 



(B.14) 



(B.15) 



As discussed in Remark 



B.l 



the function /p,n(^) behaves asymptotically as z^~p for 
large z, i.e., there exists < Cp^n < 00 such that for any £ > 0, there exists 2^0 > 
such that for any z > zq [cf. Eq. (B.16)] 



{Cp,n - e)z^-P < fp^niz) < {Cp^n + e)z^-P, 
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from which one immediately deduces 

As a consequence, for any i? < 1, one can always choose e such that the right-hand 
side is strictly smaller than 1 so that the ratio of Lp norms is then strictly positive. 



Moreover, the limiting value is 1 — R^ ^ that completes the proof of Eq. (2.9) for 
1 <p<3. D 



Remark B.l. The function fp^n{z) defined by Eq. (B.15) asymptotically behaves 
as z^~P for large z, i.e., the limit 

cp,„ = lim %M (B.16) 

exists, is finite and strictly positive: < Cp^n < oo. As for Remark [A. 1[ a rigorous 
proof of this result is beyond the scope of the paper. The asymptotic behavior from 



Eq. (2.9) for focusing modes is illustrated on Fig. A.l 



Appendix C. Analysis of Mathieu functions. 

Many algorithms have been proposed for a numerical computation of Mathieu 
functions f39j [40l [HJ |42]. The main difficulty is the computation of Mathieu char- 
acteristic numbers (MCNs). Alhargan introduced a complete method for calculating 
the MCNs for Mathieu functions of integer orders by using recurrence relations for 
MCNs [39^ . His algorithm is a good compromise between complexity, accuracy, speed 
and ease of use. Nevertheless, for illustrative purposes of the paper, we used a simpler 
approach by Zhang et al. [43 . In this approach, the problem of calculating expansion 
coefficients of Mathieu functions is reduced to an eigenproblem for sparse tridiago- 
nal matrices. We have rebuilt the computation of Mathieu functions and modified 
Mathieu functions and checked the accuracy of the numerical algorithm by comparing 
their values to whose published in the literature [131 HSl HI] • We also checked that the 
truncation of the underlying tridiagonal matrices to the size Kmax = 200 was enough 
for getting very accurate results, at least for the examples presented in the paper. 

Appendix D. Asymptotic behavior of Mathieu functions for large q. 

The large q asymptotic expansions of cen(z,g) and seri+i{z •, o) foi" ^ ^ [^^ f ) ^^^ 
n = 0,l,2, ... are [29, 45 

ce.(z, q) = Cn{q) U^^^'^^hUz) £ ^^ + e-2^^^-^/.-(z) £ ^ J (D.l) 

\ /c=0 ^ /c=0 ^ / 



k=0 ^ k=0 ^ 



where 



(cos z) 



n+1 

2n+l 



Ki--) = r^i^^^^^ii^ 



(cos Z) 



(l + sinz)^+i' 


(D.3) 


(l + sin^;)^ 
(1 -sin^)^+i' 


(D.4) 
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0.1 0.2 0.3 0.4 0.5 0, 



Fig. D.l. The functions cei(z,q) and sei(z,q) (solid and dashed lines), computed by our 
algorithm with Km ax = 200; and their approximations (circles and triangles) by the asymptotic 
expansions \DA\\D^ truncated to two terms (k = 0, \), with q = 20. 



and the coefficients Cn{q) and Sn-\-i{q) are given explicitly in [29 . The coefficients 
f^{z) can be computed through the recursive formulas given in [45j, e.g. 



f^{z) = 1, /±(z) 



2n + 1 ^ (n^ + n + l) sin z 
8cos(z)2 • 



When q is large enough, one can truncate the asymptotic expansions (D.l, D.2) by 
keeping only two terms (/c = 0, 1) and get accurate approximations for ce^ and se^+i, 
as illustrated on Fig. |D.1| 

It is convenient to define the functions 



G±(z,q) = /i+(z)±e-4^-"(^)/i-(z) + /.+ (z)^47?±e"'^''"^'^^nWE^'^^^ 



OO ^ + 

E 



^k/2 -^ "^ny-J,^ u/2 

^ k=l ^ 



in order to write 

In what follows, we will estimate the functions G^{z, q) by their leading terms, given 
that the remaining part is getting small for large q. 

Lemma D.l. For 7 G (O, f ); there exists N^ > such that for q > N. 

ZG(0,7) 



g4^W 



fe=i 



q' 



k/2 



< 



(D.5) 



Now, we establish the upper and lower bounds for the functions G^. 
Lemma D.2. Let a e (O, |), 7 G (a, |). Then, there exists N^ > such that 
for any (3 G (a, 7) and q > N^: 

\Gi{zuq)\ < ^ (l + /i;;(a)e-4^^-(-^)) V ^i G (0,a), (D.6) 

K{z2,q)\>lht{i) V^2G(/3,7)- (D.7) 
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Proof. From Lemma D.l , there exists A^^ > such that for q > N^ and zi G (0, a), 
one has 



G: 



:{zuq)\ < 2 {ht{zi) + h-{z,)e-'^^^'^^^^^) < - {l + h-{a)e 



-4:^sin{zi) 



For q > Nj and Z2 G (/3,7), one has 

K{z2,q)\ > \ {ht{z2) + h-{z2)e-'^^^'"^^^^) > i/i+(7) > 



and 



<l{ht{z2)+h-{z2)e-*^^^'"^^^^). 



G-{z2,q) - {h+{z2) - h-{z2)e-'^^^'-(^^^ 
The last inequahty imphes 

|G-(^2,9)| > min I (/i+(^2) - /i;;(^2)e-^^^"(^^)) , ^ (/i+(z2) - 3/i;;(^2)e-^^^'"(^^) 

Since h~{z2) > and h^{z2) is a decreasing function, one gets 



|Gn(^2,g)| > 2^n(7), 



that completes the proof. D 



Now we can prove Theorem 3.1 



Proof. We first consider the case i = 1. Using the symmetric properties of 
Mathieu functions ^43^, one has 



\L, ,j|P J\cen{zi,qnki)\^dzi 



J |cen(2;2,gnfei)|^a^2 



Choosing /3 = | + ^ and 7 = ^ + f , one gets 



J\cen{zi,qnki)\^dzi J\cen{zi,qnki)\^dzi 



< 



^/2 ^7 

From Lemma D.2[ there exists N^ > such that for q > N^ , 

a a 

cen{z,,q)\Pdz, = {Cn{q)r J e^P^^''^''\G+{z^,q)\Pdz, 



< {Cu{q)Y 



E 

\k=Q 



P\ \ 2^s\nziw-2k 



{h-{a)f dz^ 



< 
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where the terms e^^^^^^i were bounded by e^^^^"" for m > 0, and by 1 for m < 
(here [x] denotes the integer part of x). In addition, 



from which 

ll^nfelir 



Lp{Q\no.) 



< 



3Pa 



'b/2] 



-2py/qnki (sin /3 — sin a) l 



fe=l 



s ( 



P\( 2^sina\-2k 



{h-{a)r 



-2p^/qsin 



fc=[p/2] + l 



(Kia))' 



Taking q large enough, one can make the terms in large brackets smaller than any 
prescribed threshold e. For e = 1, one can simplify the estimate as 



\'^nkl\\Lp{n\Q^) 



<2 



3Pa 



ii^n.iiii^(^^) {7-pm{7)y 



exp 



-2p^yqnki{sinl3 - sin a) 



Substituting (3 = 7r/4-\-a/2 and 7 = 37r/8+Qf/4, one gets Eq. (3.7) after trigonometric 
simplifications. 

For i = 2, one can use similar estimates for scn+i- □ 

Appendix E. No localization in rectangle-like domains. 

Theorem |4 . 1 1 relies on the following simple estimate. 

Lemma E.l. For < a < b and any positive integer m, one has 





1 I sin(mx)|(ix > / ^iv?{mx)dx > e(a, b) > 0, 



cos{rrix)\dx > / cos^{rrix)dx > e{a^b) > 0, 



(E.l) 



where 
e(a, b) = min < 



b — a b — a 1 

4 '^ 2 



sin(n(6 — a)) 



: n = 1,2,..., 



b — a 



>0, (E.2) 



It is important to stress that the lower bound e(a, b) does not depend on m. The 
proof of this lemma is elementary. 

The proof of Theore m |4.1| is a simple consequence. 

Proof. The condition (4.2) ensures that all the eigenvalues are simple so that each 
eigenfunction is 



'^ni,...,nd \^li "") ^d) 



f sin(7rniXi/-^i)... sm{7rndXd/id) (Dirichlet), 
I cos(7rniXi/-^i)... cos{7rndXd/id) (Neumann). 
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For any open subset V, there exists a ball included in V and thus there exists a 
rectangle-like domain Qy = [<^i7^i] x ••• x [o^d^bd] C F, with < ai < bi < ii for all 
i = 1, ..., d. The Li-norm of u in V can be estimated as 



d J 

||^ni,...,nJ|Li(y) > ||^ni,...,nJ|Li(nv) =W j dXi 



|sin(7rniXi/^i)| 



CO^{TTniXi/li)\ 
''~-^ai ^ 

• -, ^ II cos(niXi)\ 



^i.-.-^d TT /" . J|sin(niXi)| ^ ^i...^, 



^ ~_i J cosfniXj TT^ t^-*- 



Tra 



% I '^% 



where the last inequality results from ([EJj). To complete the proof, one uses the 
Jensen's inequality for L^-norms and /id(V^) > jJ^di^v) = {^i — Cii)---{bd — dd) 

||'^ni,...,nJ|Lp(y) ^ ll^ni,...,nJ|Li(y)(Md(V))^"^ 1 ^ ( h - ai\^~^ ai bi 

Ti n > r~ ^^11 — 7^ ci7r^-,7r— ) > U. 

\\Unu...,nALAn) ||U„„...,„J|L^(0)(/Xrf("))^ ^ tl V ^. / ^i ^. 

Since the right-hand side is strictly positive and independent of ni, ... , n^, the 
infimum of the left-hand side over all eigenfunctions is strictly positive. D 
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